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ABSTRACT 

We study the Abelian projected £77(2) lattice gauge theory after 
gauge fixing to the maximally Abelian gauge (MAG). In order to check 
the universality of the Abelian dominance we employ the tadpole im- 
proved tree level (TI) action. We show that the density of monopoles in 
the largest cluster (the IR component) is finite in the continuum limit 
which is approximated already at relatively large lattice spacing. The 
value itself is smaller than in the case of Wilson action. We present 
results for the ratio of the Abelian to non-Abelian string tension for 
both Wilson and TI actions for a number of lattice spacings in the range 
0.06 fm < a < 0.35 fm. These results show that the ratio is between 0.9 
and 0.95 for all considered values of lattice couplings and both actions. 
We compare the properties of the monopole clusters in two gauges - in 
MAG and in the Laplacian Abelian gauge (LAG). Whereas in MAG the 
infrared component of the monopole density shows a good convergence 
to the continuum limit, we find that in LAG it is even not clear whether 
a finite limit exists. 

1 Introduction 



The dual superconductor scenario of confinement has received support from many obser- 
vations made as well in gluodynamics P3I2] as in full lattice QCD [3]. The most intensively 
investigated case was SU (2) gluodynamics. The scaling properties of many gauge depen- 
dent observables such as the Abelian string tension, the effective monopole action, the 
monopole density etc. have been checked. It has been shown that some properties of the 
monopoles in MAG can be explained by percolation theory or by free particle field theory 
[3]. Despite this progress there is lack of universality checks, i.e. the independence of 
the choice of action has not been confirmed. Apart from papers |S| and always the 
Wilson action for the gauge field has been employed. In comparison with Ref. [Sj, where 
the same TI action was used, we have much better statistics and better gauge fixing, i.e. 
lower effects of gauge fixing ambiguities. In where a different improved action was 
considered, the study was made for one value of lattice spacing only and thus no scaling 
studies were attempted. In the present paper we are aiming to make a contribution to 
closing this gap. The other, perhaps even more important problem is the dependence 
on the gauge condition used for Abelian projection. There are various opinions on this 



problem. Some authors believe that the occurrence of monopole condensation itself and, 
correspondingly, the dual superconductor properties of the vacuum have to exist and do 
exist in any Abelian gauge [7j. On the other hand recent results jH] obtained with the 
Frohlich-Marchetti monopole creation operator show that the monopole condensate de- 
pends on the choice of the Abelian gauge. In it was argued that in the Abelian gauge 
defined by diagonalization of the Polyakov loop operator the condensation of monopoles 
does not necessarily lead to formation of the Abelian flux tube between static quarks. 
We share the opinion, that any Abelian projection is made with the goal of separating 
degrees of freedom responsible for infrared physics, which thus should carry all low mo- 
menta of the original non-Abelian gauge field, from ultraviolet degrees of freedom which 
are responsible for short distance physics (if they have a sensible continuum limit) or 
might even be mere lattice artefacts. Such separation does not need to be accomplished 
in any conceivable Abelian gauge. Rather we expect that there might exist a class of 
gauges which indeed have this property. The maximally Abelian gauge MAG is a very 
likely candidate to belong to this class, and the Laplacian Abelian gauge (LAG) jTU] is 
widely considered as another good candidate. Since the analytical study of these gauges 
in the nonperturbative regime is very difficult and has not been accomplished so far de- 
spite many recent attempts (for MAG studies see, e.g. JT] and references therein), the 
numerical study is the only way to approach this problem in practice. Therefore, in this 
study we also compare some of the properties of these two gauges. From the above point 
of view, the issue of universality, it turns out that MAG is really unique to allow the 
separation of scales attempted by Abelian projection. 

The paper is organized as follows. In section |21 we specify the technical tools, in 
particular the improved action used here in contrast to the Wilson action and the method 
of gauge fixing. Then, in section El we briefly report on the evaluation of the string 
tension for the purpose of calibrating the lattice scale corresponding to the improved 
action. Section 14.11 contains our observations concerning the scaling properties of the 
monopole densities and their IR and UV components for both actions in the maximally 
Abelian gauges. In section 14.21 we show that similarly to MAG the monopole clusters 
obtained in the Laplacian Abelian gauge might be splitted into IR and UV components 
but their scaling properties are quite different from those observed in MAG. Section is 
devoted to the Abelian dominance study. Our results indicate universality of the Abelian 
dominance in the continuum limit. Finally we summarize our findings in section |3 

2 Simulation details 

To address again the question of universality, we employ here the tree level improved 
action of the form [T2] 



where S p i and S rt denote plaquette and 1x2 rectangular loop terms in the action, 



the parameter uq is the input tadpole improvement factor taken here equal to the fourth 
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In our simulations we have not included one-loop corrections to the coefficients, for 
the sake of simplicity and also to be able to compare with the results of Ref. [3] after 
making a few improvements in comparison with this work in other directions. First, we 
improved substantially the gauge fixing as will be explained later. Second, we have now 
better statistics and worked on larger physical volumes. This has allowed to determine 
more reliably various Abelian observables and their infrared part. Third, we have used 
a new smearing techniques which enabled us to make more precise measurements of the 
non- Abelian string tension. This was necessary to assess the nonperturbative scaling of 
various monopole densities. 

We also make a comparison of these Abelian observables obtained in MAG and LAG, 
respectively. The MAG is fixed by the maximization of the lattice functional 

F W = E Tr M^st^U) , (3) 

with respect to local gauge transformations 

For MAG we applied the simulated annealing algorithm. The details of the gauge fixing 
procedure can be found in We have applied the algorithm to 10 randomly replicated 
gauge copies of each Monte Carlo configuration in the hope to find among the 10 local 
maxima one closest to the global maximum. This procedure proved to be the best so far 
to fix MAG as well as to fix center gauges [Tl] . Although there is no proof we hope that 
our results for gauge noninvariant observables are numerically close to those we would 
obtain evaluating it at the gauge equivalent representant carrying the global maximum of 
(J3J) for every configuration of the gauge field. 

Fixing the LAG amounts to finding the eigenvector with the lowest eigenvalue of the 
covariant Laplacian operator in the adjoint representation, 

D nm — ^ (2S nm 5 i? n ^5 mjn _|_^ R n _£ L fl 8 mtTl ~{i) (5) 
with the adjoint link variable 

= Tr (a a U n>fi a b UlJ . (6) 

The gauge transformation g n is then determined by the requirement to rotate this eigen- 
vector 0^ to the 3 rd color axis at every site n: 

3 , — 

Pn&3 = 0« 9n(T a gL Pn = \/ <\>\ • ( 7 ) 

0=1 

The simulations with the action JI} have been performed with parameters given in 
Table ^ The parameter uq has been iterated over a series of Monte Carlo runs in order 
to match the fourth root of the average plaquette P. The corresponding entries give an 
impression of the accuracy of matching. For two values of f3i mp we simulated lattices of 
two sizes. The value of the parameter uq fixed on smaller lattice was used as an input for 
larger lattice simulations. The string tensions obtained on these lattices are in agreement 
within error bars. The smaller (larger) lattices were used to study LAG (MAG). 



Table 1: Details of the simulations with improved action 
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Table 2: Details of the simulations with Wilson action 
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3 The non-Abelian string tension 

In order to fix the physical lattice scale we need to compute one physical dimensionful 
observable the value of which is known. For this purpose we choose the string tension 
a. The string tension for action (0) was computed long ago in J5 but we will improve 
this measurement according to present standards. We use the hypercubic blocking (HYP) 
invented by the authors of Ref. [H| to reduce the statistical errors. This method has been 
successfully applied to static potential calculations in SU(3) gluodynamics [13 HH1 UZj 
and in full lattice QCD at finite temperature |18j . After one step of HYP, about 20 
sweeps of APE smearing [HI! were applied to the space like links. The spatial smearing is 
made, as usually, in order to variationally improve the overlap with a mesonic flux tube 
state. In Fig. ^ we compare potentials obtained with and without HYP procedure. As 
was observed in the cited above papers the HYP potential differs essentially by a constant 
shift corresponding to reducing the static source self-energy. One can see from the figure 
that HYP decreases both statistical errors and effects of rotational invariance breaking. 
Since HYP changes the potential at small distances we included only distances R/a > 2 
into our fits of the static potential. The resulting values for Vera 2 are also included in 
Table H The string tension was also calculated with Wilson action. In this case APE 
smearing for space links and the additional trick of link integration [20] for time links were 
used in the evaluation of Wilson loops. The results for v aa? and details of simulations 
with Wilson action are presented in Table El 




Figure 1: Non-Abelian potentials for the TI action obtained without (circles, SM) and 
after (squares, HYP) hypercubic blocking (both with spatial smearing) vs. R/a for T/a = 
5 at f3i mp = 3.4. 

4 The monopole density 
4.1 MAG 

After fixing the Abelian gauge the Abelian projection can be made: 

Un,n Cn,/i M n,(i j (8) 

where the Abelian field is contained in u n ^ = diag(e* en ' M , e~ t6 "^), 9 n ^ G (— n, n], and C n ^ 
is the coset field describing charged gluons. The Abelian plaquette angle 

@n,fii/ &ifin,v dv@n,fM (9) 

is decomposed into regular and singular parts: 

9 n ,iiv = On^u + 2irm n ^ , 9 n ^ v E (-7T, 7r] (10) 

@n,fj,v is a physical Abelian flux through the lattice plaquette {n, /xz/}, and m ntfiu counts the 
number of Dirac strings through this plaquette. The magnetic currents are then defined 
as follows: 

k n ,[i = -^iivafidydn^p = — — t fJiUa j3d u m nta i3 (11) 
We will measure the monopole density in lattice units p lat defined as 
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Since monopoles are three dimensional objects their physical density is related to the 
lattice density by a 3 p = p lat . With Wilson action, the first measurement of the monopole 
density in MAG gauge has been made in j2Ij with participation of the authors of the 
present paper. That result was interpreted in the sense of asymptotic scaling. In fact, 
the observation of asymptotic scaling at (3 values in the range from 2.4 to 2.6 would seem 
rather strange today. It is known that the string tension does not follow the two-loop 
renormalization group formula over this range in (3, such that the result obtained in [21] 
actually implies the divergence of the monopole density with f3 — > oo. The situation has 
been partially clarified by Hart and Teper [22] |2Hj . These authors found that on large 
enough lattices the network of magnetic currents in each configuration consists of one large 
cluster and many other clusters with much smaller size. The spectrum of cluster sizes falls 
into two very distinctive parts, disconnected by a gap. The density of currents forming 
the largest (percolating) cluster p IR has been measured in units of the string tension 
and a first indication of scaling of the ratio pm/cr 3 ^ 2 has been found. More accurate 
measurements [21] have corroborated this kind of scaling behavior. The continuum limit 
for this ratio was determined as pm/cr 3 ^ 2 = 0.65(2). 

Another important result obtained in Refs. [221 122] was the observation that the 
largest cluster alone produces almost the full monopole string tension. This fact has 
allowed the authors of Ref. [23] to call the monopoles belonging to this cluster "infrared 
monopoles" (IR) while the monopoles from the remaining clusters were called ultraviolet 
(UV) monopoles, implying that these monopoles are not relevant for IR physics. We 
should mention that, despite the fact that they are not relevant for the confining proper- 
ties of the vacuum as supported by numerical observation, their relevance for the topology 
and therefore for the chiral properties of the vacuum has not yet been explored. We will 
keep (and have already used) the above notation, quoting IR and UV monopoles in the 
following. In Ref. [21] it has been demonstrated that the density of UV monopoles, hence 
the total density, diverges in the continuum limit. 

It should be noticed that in [221 I2H] a single, supposedly percolating, cluster with 
a size much larger than all other clusters in the given configuration was only found on 
large enough lattices. For decreasing lattice size L and fixed lattice spacing this gap 
disappears, i. e. the largest and the second-largest cluster are of similar size. Consequently, 
the important property that the largest cluster alone produces almost the full monopole- 
related string tension, is lost. This implies "splitting" of the largest cluster when the 
lattice volume decreases. It was also found in f2"2~] I2H] that the critical value of the lattice 
size L crit , below which the largest cluster splits, is in the range of lattice sizes, where 
physical quantities do not show large finite volume effects, and, moreover, L crit , measured 
in physical units, increases with decreasing lattice spacing. Such behavior implies that in 
the continuum limit the gap in the spectrum of cluster sizes might disappear and a clear 
separation of IR and UV monopoles may become impossible. 

The solution of this problem was suggested in [21] . It was found empirically that the 
splitting of the largest cluster leads to formation of clusters (rare for large enough lattices) 
with nonzero winding 



Such clusters might be very large on given configurations or be of moderate size. In both 
cases they extend through the whole lattice at least in one direction and thus should 
be considered as relevant for the infrared physics. Furthermore, two or more of such 




(13) 



clusters (forming together a combination of clusters with zero total winding w M ) might 
form boundaries of the same Dirac sheet, which is closed in one or more directions due to 
periodic boundary conditions. Let us note that in case of two wrapping clusters present 
in one configuration they unambiguously form the boundary of the same Dirac sheet. 
When three wrapping clusters are present they also form the boundary of one Dirac 
sheet or boundaries of two Dirac sheets. In the latter case one of the clusters forms part 
of the boundaries of both Dirac sheets while two others form part of the boundary of 
one of those Dirac sheets. It is natural to consider such clusters as one cluster when it 
concerns the determination of clusters relevant for IR physics. It is also clear that the 
splitting phenomena can be, at least partially, ascribed to the annihilation of parts of the 
percolating cluster through boundary conditions leading to formation of two disconnected 
clusters which still form a boundary of one Dirac sheet. 

Based on these observations, it was suggested in j21j to define, for each configuration, a 
single IR cluster as the union of the largest cluster (which might have trivial or nontrivial 
winding) and all clusters with nonzero winding. Numerical evidence was further presented 
in |21j that under such definition the size of the largest cluster changes smoothly with 
the lattice size for physically large lattices. There are preliminary results of ours to 
be published elsewhere, showing that with this definition the IR cluster alone reproduces 
almost the full monopole-related string tension. We will use this definition in what follows. 

As for lattice fields generated with TI action (|Tj) . the monopole density has been 
measured by Poulis 5j. He concluded that the total density has correct scaling in the 
continuum limit. In the light of the discussion above this would mean that the TI action 
would take an exclusive role. However, as we show below, this conclusion was wrong. 
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Figure 2: The monopole cluster length distribution N(l) at (5 = 2.55 for the Wilson action 
(left) and at (3 imp = 3.5 for the TI action (right). 

In Fig. El we show the distributions of the monopole clusters length for Wilson action 
at (3 — 2.55 and for TI action at (3i mp = 3.5. Note, that a single IR cluster per configu- 
ration has entered the histograms which has been defined for the configuration at hand 
as described above. One can see that these distributions are qualitatively similar, i.e. for 
the TI action we also observe a clear splitting of the clusters into IR and UV clusters. 

Our results for the densities of IR and UV monopoles and the total density, taken in 
lattice units, are presented in Tabled In physical units of cr 3 / 2 , the results are shown 
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Figure 3: The monopole densities in MAG for the case of the TI action. The dashed line 
shows a quadratic fit. 

as a function of lattice spacing (in units of V ad 2 ) in Fig El One can see that the IR 
density converges to a finite value in the limit a —>■ 0. In contrast to this, the density 



Table 3: Density of monopoles 
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of UV monopoles, and thus the total density, behaves divergent in the continuum limit. 
We thus find that the TI action leads qualitatively to the same picture as was observed 
before with the Wilson action. 

To make a quantitative comparison we plotted in Fig. |U p IR (left) and puy (right) 
for two actions. From these figures one can see quantitative differences. The continuum 
value of pm/cr 3 / 2 , obtained with a quadratic fit, is 0.50(1) for TI action and 0.71(2) 
for Wilson action, i.e. they differ by factor 1.4. We should make a comment on the 
different procedures of calculation of Pir/u 3 ^ 2 for Wilson action in the present paper and 
in Refs. |24( . One difference is that in these papers only subsets of the full ensembles 
of gauge field configurations used in the present paper were employed. Another, more 



important difference is that, in these earlier papers values of a from literature were used, 
while here we are using values of a calculated, as described in the previous section, i.e. 
on the same set of configurations on which the monopole density was calculated. 

The observed difference in pm measured for TI and Wilson actions, though not being 
drastically large, means that the present definition of IR density is not universal. This 
makes it difficult to ascribe to it a meaning as a physical, gauge invariant density. It is 
evident that the source of the discrepancy in the values of pm might be the appending of 
UV monopoles, i.e. small loops, to IR monopole clusters. Since the TI action suppresses 
UV degrees of freedom stronger than Wilson action, it is natural to expect that this 
additional length assigned to the IR cluster is smaller for TI action. Whether this is the 
only reason deserves further investigation. 

The density of UV monopoles is reduced much more substantially, roughly by a factor 
2.5, as can be seen from Fig. |U (right). We can say that the TI action indeed suppresses 
(part of) the UV degrees of freedom. As generally expected for an improved action, one 
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Figure 4: Comparison of the monopole densities obtained with TI (full symbols) and 
Wilson actions (open symbols) in MAG: left - IR monopole density; right - UV monopole 
density. 

can also see earlier and faster convergence to the continuum limit. 

As it has been mentioned above, the UV monopole density diverges. It is natural 
to ask to which power of a -1 this divergence is compatible. It was first found in [23] 
that for the Wilson action puy ~ 1/a. In Ref. [2E] this was confirmed with higher 
confidence. Fig. 0] (right) shows apuyja for both actions. For TI action this ratio seems 
to rapidly converge to a finite value in the continuum limit as soon as V a a? < 0.25 . 
This implies that puv ~ 1/a also for TI action. Note, that the convergence for the TI 
action is faster than for the Wilson action, as can be seen from Fig. 0] (right). If the 
existence of a reasonable continuum limit in the latter case would be confirmed, then only 
for V 'a a? << 0.1, and the limit value would be markedly larger than for the improved 
action. In any case, the data clearly show that the UV monopole density and therefore 
the total density of monopoles is not universal. 

The TI improved action has corrections of order 0(a s a 2 ) and 0(a 4 ) while the Wilson 
action has corrections of order 0(a 2 ). Thus it is natural to expect that some contribution 
of lattice artifacts to the monopole density is suppressed in case of the improved action. 




Figure 5: The small (UV) clusters' length distribution in MAG for the Wilson action at 
various f3 (left) and for the TI action at various (3 imp (right). Curves are fits to Eq. (fTlj) 

Now we are able to conclude that a considerable part of the UV monopoles in the Wilson 
action case (more than 50 %) are lattice artifacts. In contrast to this, we notice that IR 
monopoles are not much affected by lattice artifacts 1 . Whether the value of pm obtained 
with the improved action is already the final universal one is an open question. This 
should be checked in simulations with other improved actions. 

The interesting, yet unanswered question is the physical role of UV monopoles. It was 
found in [22 I2H] and then confirmed with higher precision in that the number N(l) 
of small clusters of length I falls like 

^r = cC9)/r, (14) 

where 7 ~ 3. The value 7 = 3 was shown to be in agreement with percolation theory and 
also to be derivable within the polymer approach to the field theory of free or Coulomb- 
like interacting scalar particles jl]. Our data for both actions are also in agreement with 
relation (fUJ). with values of parameter 7 close to 3, as can be seen from Fig.Eland TablelH 



Table 4: Parameters of the fits to Eq. (fT4^) of the small clusters's length distribution in 
MAG for both actions. Respective fit ranges (l m i n ,l m ax) are also shown. 



p 


7 


c(P) 


I min 


^max 


Pimp 


7 


c([3) 


Imin 


Imax 


2.45 


2.93(2) 


0.23(2) 


10 


70 


3.3 


2.97(4) 


0.12(1) 


8 


40 


2.50 


2.99(2) 


0.20(1) 


12 


50 


3.4 


3.01(7) 


0.08(2) 


10 


30 


2.55 


3.08(2) 


0.17(1) 


10 


50 


3.5 


2.87(9) 


0.03(1) 


10 


30 


2.60 


3.11(6) 


0.12(2) 


14 


40 













Let us come back to observation that puy ~ -. This implies that magnetic currents 
from small clusters have a finite density per unit of 2D volume rather than per unit of 

1 The IR density is definitely affected by the roughness of monopole currents inside the largest cluster, 
and hence not completely free of discretization artefacts. 



3D volume which is actually the case for the IR magnetic currents. It has been recently 
verified that the density of P-vortices in the indirect Z{2) center gauge is finite in the 
continuum limit (2Z1- On the other hand, it is known that in this gauge P-vortices and 
monopoles are highly correlated |28*1 It is then natural to assume that magnetic 
currents belonging to the small clusters "populate" the surfaces formed by P-vortices 
with some constant density. The strong reduction of the density of UV monopoles in the 
case of TI action in comparison with Wilson action suggests that the density of P-vortices 
will be suppressed, too. This should be checked in a future calculation. 

We now introduce a regularized UV monopole density by summation of the average 
number of small clusters of length I per lattice volume, i.e. N(l), multiplied by the length 
I. In this definition we shall exclude clusters below a certain (minimal) length scale l p h 
specified in natural units: 

Pvv{l P u) = ^3 £ N(l)l, (15) 

l>l p h/a 

We call this UV monopole density "constrained density". This definition counts all 
monopole currents in small clusters with a length above or equal l p h, and the emerging 
density depends on it as a parameter. Thus, very small loops sensitive to the ultraviolate 
cutoff are excluded. In Fig. |H]we show puv(lph) /cr 3 ^ 2 as a function of l p h for both actions. 
One can see that scaling is very good as long as l p h\f® < 2 in the case of the TI action 
and l v h\fo < 4 in the case of Wilson action. In general, scaling becomes worse at larger 
values of l p h\fo '. This might be the consequence of the IR cluster splitting discussed above 
which underlies the splitting into UV and IR monopoles: some large clusters which actu- 
ally should belong (are akin) to IR monopole clusters were counted as small ones because 
of trivial winding. Although such clusters are relatively seldom their number increases 
with increasing (3. To exclude the effect of these ambiguously identified "UV" clusters we 
subtract the contribution of clusters with I > l p h = c/ y/a and plot in Fig. [7| the difference 

Puv{l P h) - Puv(Jph) (16) 

with coefficients c ~ 13.6 and c ~ 10.65 for Wilson action and TI action, respectively. 
One can see that now scaling is uniformly well satified for all lower cutoffs l p h < l p h except 
very small ones. Thus we come to the conclusion that the UV monopole density derived 
from the small cluster density N(l) as defined in eq. I|15|L i.e. when clusters close to the 
cutoff scale are excluded, shows good scaling, i.e. is independent of a similar to the IR 
monopoles density (derived from the IR clusters). 
Using eq. (fT3)l with 7 = 3 we get 

Puv(l ph ) = ^'(l ph /a), m = j£j (17) 
We fitted the data in Fig. [7| by the function 

ciV>'(W°) + c 2 (18) 

and found that this function describes well our data, especially for the TI action case, 
with constants q only weakly dependent on f3. This implies that for small a (large l p h/o) 
the density behaves as puv(lph) ~ &/lph since ip'(z) ~ 1/z for large \z\. This fact is 
in agreement with 1/a behavior of puv we argued for, since it can be unambiguously 
regulated by the assignment puv(l p h = 4a) — > puv- 




Figure 6: The "constrained" UV monopole density derived from the small clusters' length 
distribution Eq. (JTKJ) for TI (left) and Wilson (right) actions. 




Figure 7: The "constrained" UV monopole density derived from the small clusters' length 
distribution and corrected for ambiguous clusters according to Eq. (fTEJ). for TI (left) and 
Wilson (right) actions. The curves show fits by Eq. (jTSJ). The data were lifted by a 
constant for better readability of the figures. 



4.2 The Laplacian Abelian gauge (LAG) 

To calculate the lowest eigenvector of the covariant adjoint Laplacian operator (J3J) we used 
the Arnoldi algorithm [3U|. This algorithm, as well as others used to solve this problem, 
requires large memory increasing fast with lattice size. Therefore, the measurements in 
LAG for Wilson action have been made on smaller lattices than shown in Table 121 for 
(3 = 2.45, 2.55, 2.6, but with large enough physical size, which was never smaller than 1.4 
fm. In case of the TI action we made measurements in LAG only for four values of the 
coupling constant, because this proved to be enough for our purposes of comparison with 
MAG. 

We first present the cluster length distribution for the two actions in Fig. |H1 One can 
see that in LAG the separation into IR and UV clusters works very well. We note that 
there are only rare cases of clusters with nontrivial winding even for most fine lattices. 
Our results for various densities are presented in Fig. El We found that the total monopole 
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Figure 8: The monopole cluster length distribution N(l) in LAG at (3 = 2.6 for the Wilson 
action (left) and at (3i mp = 3.5 for the TI action (right). 

density in LAG is substantially higher than in MAG in agreement with earlier observations 
[TO] made for Wilson action. We further looked at IR and UV densities separately and 
found that increasing of the density is true for both of them. For our finest lattice 
these densities in LAG are 2-3 times higher than in MAG. It is not clear from our data 
whether the IR density in LAG converges to a finite value as it is the case in MAG. 
This is most probably due to a substantial (more fractal) admixture of lattice artefacts 
to the IR clusters. We also observed that the UV monopole density puy is much more 
strongly diverging in LAG than in MAG and is compatible with puy ~ I /a 2 rather than 
with puy ~ 1/a, as observed in MAG. Accordingly, the constrained density puvi}ph) was 
found to be divergent as The small clusters' length distribution was fitted to eq- lpPfj) 
with a resulting parameter 7 in the range 2.7 — 2.8, i.e. significantly smaller than for 
MAG. 
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Figure 9: Various monopole densities in LAG for TI action (left) and for Wilson action 
(right). Lines are drawn to guide the eye. 



5 The Abelian string tension in MAG 

To estimate the Abelian string tension we calculate the Abelian potential V ab (R) using 
(spatially) smeared Abelian Wilson loops W ab (R,T). As usually V ab (R) is defined as a 
limit 

V ab (R)=\im Vab(R,T), (19) 

T — >oo 

where the potential estimator V a b(R, T) is 

t/ It) t\ i W ab {R,T + a) 
aV ab (R,T) = -\o g WMT) . (20) 

It is important to check that V ab (R) is unique, i.e. independent of the operators used to 
create the Abelian flux tube. One can get different such operators varying the number 
of smearing sweeps N sm . In Fig. [TUlwe show V ab (R,T) for N sm = 3 and 100. For large 
number of sweeps the behavior of V ab (R,T) clearly shows absence of the positivity. For 
small number of sweeps the behavior is similar to the case when positivity is fulfilled. This 
is presumably due to higher excitations: in this case they are not suppressed for small 
T/a. Lack of positivity for gauge dependent correlators in covariant gauges was discussed 
recently in [31]. The most important for us observation to be read off from Fig. ITU1 is 
that at large enough T results agree with each other. This implies that V ab (R) is indeed 
defined uniquely. 

We found that the behaviour of the ratio cr ab /a for TI action is qualitatively similar to 
that for the Wilson action, see Fig. El For both actions the ratio is between 0.9 and 0.95 
for all considered values of lattice spacing. We thus definitely confirm the universality of 
Abelian dominance in the continuum limit. Our results, due to large statistical errors, 
coming mainly from the determination of the non- Abelian string tension, do not allow to 
determine precisely the continuum limit of the ratio u ab /a. 
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Figure 10: The Abelian potential V a b(R,T) in MAG at f3i rnp = 3.3 vs. T/a for R/a = 4 
for 3 (open circles) and 100 (full circles) smearing sweeps. 
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Figure 11: Ratio between the Abelian string tension (in MAG) and the non- Abelian string 
tension vs. lattice spacing for TI action (left) and for Wilson action (right). 



6 Summary 



In this paper two important questions on the properties of the Abelian projection were 
addressed: universality and gauge dependence. Comparing results obtained with Wilson 
and TI actions on lattices with varying lattice spacing we confirmed that the Abelian 
dominance passes the universality check in MAG. Moreover, this universality holds in the 
continuum limit. We found that in the continuum limit the ratio <J a b/<J seems to be in 
the range from 0.90 to 0.95. No convergence to 1 was observed contrary to our earlier 
results |21] seen with smaller statistics. We have not yet accomplished the check of the 
monopole dominance universality for the string tension although our preliminary results 
confirm it. They will be published elsewhere. 

For the monopole density we found qualitative similarities: for both actions pm is 
finite in the continuum limit, while puv is divergent as l/a. On the quantitative level we 
found a violation of universality for the densities. This implies that an UV contribution 
might be substantial in the measured pm, or in other words, pm needs to be properly 
renormalized. 

We further introduced a constrained UV density puv(lph), ec L- ((HI) which is determined 
by counting only monopole loops longer than some physical length l p h. We found that 
this density scales properly in the continuum limit: 

CO~ 

Puv(l P h) = i— , (21) 

where the coefficient c is independent of a but has a different value for Wilson and TI 
actions, i.e. is non-universal. 

A comparison of the monopole densities in MAG and LAG, made for both actions, 
revealed that both pm and puv are 2-3 times higher in LAG. It is not clear from our 
data whether pm in LAG is finite in the continuum limit. The UV component p\jy 
behaves like ~ l/a 2 in LAG contrary to the divergence ~ l/a found in MAG. This adds 
to other doubts expressed in the literature about the usefulness of the LAG. Thus, 
the maximally Abelian gauge turns out to be particularly suited to the separation of IR 
and UV degrees of freedom. 
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